# Presented By ZhaoEnxing
# time: 2021/10/20 10:10
# 根据长期实验,飞机的最大飞行速度服从正态分布。现对呆新型飞机进行了15次试飞,测得各次试飞时的最大飞行速度(米/秒)为:
# 422.2 417.2 425.6 425.8 423.1
# 418.7 428.2 438.3 434.0 412.3
# 431.5 413.5 441.3 423.0 420.3
# 试对该飞机最大飞行速度的数学期望值进行区间估计(置信概率0.95)。
"""
总体服从正太分布，小样本，n=15 sigem不知道
"""
import pandas as pd
import numpy as np
from scipy import stats
n = 15
dataSer = pd.Series([422.2, 417.2, 425.6, 425.8, 423.1,
                     418.7, 428.2, 438.3, 434.0, 412.3,
                     431.5, 413.5, 441.3, 423.0, 420.3])
sample_mean = dataSer.mean()  # 计算平均值
sample_std = dataSer.std()  # 计算标准差
print('样本平均值=', sample_mean, 'm/s')
print('样本标准差=', sample_std, 'm/s')
df = n - 1
alpha = 0.05
print("显著性水平:{}".format(alpha))
tci = stats.t(df).isf(alpha / 2)
print("对应t值：{}".format(tci))
se = stats.sem(dataSer)
print(se)
a = sample_mean - tci * se
b = sample_mean + tci * se
print('在95的置信水平下，飞机最大飞行速度的数学期望值进行区间估计 CI = (%f,%f)' % (a, b))
